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^ ■ Abstract 



We study the absorption of a minimally coupled scalar in the gravitational background 
created by a stack of near-extremal black three-branes, and more generally by M2, M5 and Dp 
branes. The absorption probability has the form P(l) = P (l)fi(X), where Pq(1) is the partial 
wave's absorption probability in the extremal case, and the thermal factor /;(A) depends on the 
ratio of the frequency of the incoming wave and the Hawking temperature, A = lo/ttT. Using 
Langer-Olver's method, we obtain a low-temperature (A ^> 1) asymptotic expansion for P{1) 
with coefficients determined recursively. This expansion, which turns out to be a fairly good 
approximation even for A ~ 1, accounts for all power-like finite-temperature corrections to Pq{1), 
and we calculate a few terms explicitly. We also show that at low temperature the absorption 
probability contains exponentially suppressed terms, and attempt to develop an approximation 
scheme to calculate those. The high-temperature expansion is also considered. For the s-wave, 
the low-temperature gravity result is consistent with the free finite-temperature field theory 
calculation, while for high temperature and higher partial waves we find a disagreement. As a 
check of the approximation methods used, we apply them to the Dl — D5-brane system, and 
compare results to the known exact solution. 



1 Introduction 



It is believed that a strongly coupled gauge theory at finite temperature can be studied us- 
ing its gravity dual description in the framework of the AdS/CFT correspondence. Exam- 
ples relevant for d = 4 gauge theories include Witten's interpretation of the confinement- 
deconfinement phase transition in terms of the Hawking - Page transition in AdS and 
calculation of the free energy of the finite-temperature M = 4 SYM which matches (up 
to a factor of 4/3) the free energy of the near-extremal three-branes (for review and 
references see section 3.6 of [0]). 

Calculation of the absorption cross section by .D3-branes revealed a remarkable agree- 
ment between conformal field theory and gravity results, obtained, correspondingly, at 
weak and strong 't Hooft coupling 0. The agreement has been explained (at least for 
the s-wave absorption) by establishing a non-renormalization theorem enjoyed by certain 
J\f = 4 SYM correlation functions |§, and serves as one of the few direct checks of the 
conjectured AdS/CFT correspondence. Supersymmetry of the D3-branes' world- volume 
theory is crucial for the existence of the above mentioned non-renormalization theorem. 
At finite temperature the supersymmetry is broken, and one cannot expect the results 
computed at strong and weak coupling to agree. 

However, for the Dl/ D5 model of a five-dimensional black hole the agreement between 
greybody factors and the entropy calculated in two-dimensional CFT and in gravity 
holds even in the absence of supersymmetry. The absorption cross-section has the form 
a = aof(uj,TL,Tji), where <7o is the cross-section at zero temperature and the thermal 
factor /(a>, T L , T R ) can be computed either from (super)gravity (in the low energy regime, 
i.e. when ujt\ <C 1, ur 5 <C 1) or by considering a finite-temperature two-point function 
in the dual CFT ||, 0, ||. A perfect matching between CFT and gravity results can be 
explained by arguing that at sufficiently low energies the excitations of the D1/D5 system 
are correctly described by the moduli space approximation even at strong coupling 0. 

A priori, the existence of a non-renormalization theorem not related to supersymme- 
try in d — 4 can be considered as an extremely remote possibility, and indeed in this 
paper we do not find an exact agreement between gravity and weakly coupled field the- 
ory calculations. However, the details of the disagreement are intriguing enough, and in 
the spirit of AdS/CFT we regard our results obtained from gravity as a prediction for a 
strongly coupled finite-temperature gauge theory. 

For the Dl-D5-brane system, it is fairly straightforward to calculate the thermal factor 
f{oj, Tl, Tr) on the gravity side, since in the low energy regime the exact solution of the 
corresponding wave equation is known. This is not the case for the generic black brane 
solution, and for the three-brane solution corresponding to finite temperature d = 4 gauge 
theory in particular. It is a purpose of this paper to develop a reliable approximation 
scheme which would allow one to calculate greybody factors for the case of a non-extremal 
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three brane, and more generally for M2-, M5- and Dp-branes. Earlier attempts to 
calculate the absorption probability for non-extremal three-branes were made in ||, |9|], 
[|IIJ. However, we will show that only the zero-temperature limit of the results obtained 
there is correct. 

In the low energy regime (uR <C 1) the radial part of the wave equation in a typ- 
ical black-brane background has several regular singular points. Three of them always 
correspond, respectively, to the position of outer, inner horizons and spatial infinity, 
and the rest lies on a unit circle in the complex plane. Since the number of singular 
points exceeds three, the exact solution is unaccessible. However, the equation con- 
tains a parameter A = u/nT making it possible to study "low-temperature" (A 3> 1) 
and "high-temperature" (A <C 1) asymptotics. While obtaining the high-temperature 
asymptotics is a straightforward application of perturbation theory, the low-temperature 
expansion requires a careful WKB analysis of the differential equation in the presence of 
singular points. Basic principles of such an analysis have been developed a long time ago 
11], [l^] (see also [jn| and a classic book by Olver fll4l). Using Langer - Olver's method, 
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we are able to obtain uniform asymptotic expansions (with coefficients determined recur- 
sively) for the solutions of the wave equation. The uniformity of approximation is vital 
for the correct calculation of the absorption probability. Olver's theorem guarantees this 
property only when the coefficients and the corresponding error-control function con- 
verge everywhere on the interval of consideration. For non-extremal branes this interval 
always contains two singular points (a situation more general than the one originally 
treated in fllT , [T2||), and therefore the convergence should be checked in each case inde- 



pendently. Once the convergence has been established, all power-like corrections to the 
zero-temperature (extremal) result can in principle be determined using the recursive re- 
lations. We obtain explicit asymptotic expansions for the black three-brane background 
(Section |3.1| ) and for non-extremal Dl — D5-brane system (Section ^), and demonstrate 
(in Section |]) that all M- and black-brane backgrounds can be treated using this ap- 
proach. Being an asymptotic expansion in powers of 1/A, the method cannot account 
for the contributions of the form exp (—A) (and for the s-wave those are the only ones 
present). In Section |3.2| we make an attempt to build an approximation scheme based 
on the hypergeometric function which would allow the calculation of these exponentially 
small contributions. We find the leading term of the corresponding expansion, but the 
subleading corrections are controlled by an integral equation which we were unable to 
solve. Section |3]3| deals with the high-temperature expansion. In Section |5] we com- 
pute the absorption cross section in field theory using the method of finite-temperature 
Green's functions. Appendices |A], [FJ, ^ contain some relevant technical details. 
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2 Summary of results 



Our main result can be formulated as follows: the absorption probability for a min- 
imally coupled massless scalar in the gravitational background created by a stack of 
near-extremal black three-branes can be written as P(l) = Po(l)fi(\), where 

4tt 2 /^\ 4 '+ 8 , , 



(Z + l)! 2 (Z + 2)! 2 V 2 

is the absorption probability for the extremal D3-branes. The absorption cross section 
(greybody factor) for the Z-th partial wave is related to the absorption probability Pi via 

_ 2^T(l + d-2) (Z + d/2-1) 

^-!r(rf/2-i/2)r(/ + i) w ' 1 J 

where eZ is the number of dimensions transverse to the brane [ 15] . For the thermal factor 
fi(X) we obtain a low-temperature (A = uo/nT ^> 1) asymptotic expansion 

m = 1 + M^-i)(^-4) 

2z/(z/ 2 - l)(z/ 2 - 4)(z/ - 3)(9i/ 4 - 23z/ 3 - 114z/ 2 + 188v + 600) 
+ 225A 8 

+ O (l/A 12 ) , (2.3) 

where u = I + 2. The function //(A) for A ^> 1 also contains exponentially small terms 
(see Q3.39Q ). The high-temperature limit of //(A) for various I is given by ( ft.47| ), ( |3.48|) , 
( |3.49|) . In particular, for the s-wave the gravity calculation gives 

(1 + (e~ A ) for A > 1, 

/(Wav(A) = g / (tt 2 + 12(2 -log 2) log 2) A 2 \ r , , ( 2 - 4 ) 

+ 1 48 + -J f ° rA<<1 ' 

while in field theory we get /o,fi e idth.(A) = coth(7rA/4), i.e. 

(l + 0(e- x ) forA>l, 
fo, field th .(A) = \ 4 / ^2 A 2 \ (2.5) 

u( i+ forA<<L 

For Z = 1 the disagreement is visible already at low temperature, 

/(Wav(A) = l + 48/A 4 + 0(e" A ) , (2.6) 

/o,fiew th .(A) = l + 8/A 2 + 48/A 4 + 32/A 6 + 0(e- A ) . (2.7) 

As mentioned above, we do not expect free field theory and gravity results to agree, 
especially in the regime w<T, where the resummation of an infinite subset of diagrams 
in the finite-temperature gauge theory is required. We rather regard ( |2.4j ) as a prediction 
for a strongly coupled gauge theory at non-zero temperature. 
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3 Absorption by non-extremal black three-branes 



The metric of the black three-brane is given by |16] 

t 2 + dx ■ dx /—— 



where 

H(r) = l + ^, /(r) = l-^. 

The position of the inner horizon corresponds to r = 0, the outer horizon is located 
at r = r . We are interested in the region of the parameter space characterized by 
ujR <C 1, TR <C 1 with the ratio A = u/ttT remaining arbitrary. The condition TR <C 1 
implies that the metric we consider is near-extremal (ro -C R). For ro <C R the Hawking 
temperature is related to the non-extremality parameter ro by T = tq/ttR 2 . The horizon 
area is given by A H = n^r^ H(r ) — > -K^r^R 2 for r <C R. 

Using p = ur, p = ur Q , one can write the radial part of the Laplace equation for a 

minimally coupled massless scalar as 

d 2 <j> 5p 4 - pi d<j> _ pH(l + 4) A p\l + {uRfl _ 

dp 2 + p(p 4 - Pi) dp p^-pt^ ( P *-Pi) 2 [S - 2) 

In the extremal case ro = 0, equation ( |3.2| ) in terms of z = r 2 reduces to 

4z 3 <P" + 12zV + (u 2 R 2 - 1(1 + 4)z + u 2 z 2 ) = 0, (3.3) 

solutions of which are knownQ Exact solutions of (|3.2| ) are not known, and therefore one 
must seek a reliable approximation scheme. Introducing x = r 2 /r 2 , we can write ( |3.2| ) as 

_ (1 + x 2 ) d<p _ 1(1 + 4) | (\ 2 + p 2 /x 2 ) ^ = Q 
dx 2 x(l — x 2 ) dx 4x 2 (l — x 2 ) 4x(l — x 2 ) 2 

Following we shall calculate the absorption probability to the leading order in uR <C 1 
by introducing the "inner" and "outer" regions for the equation ( |3.4j ), and matching the 
corresponding solutions. In the outer region, p ^> p , p ^> (ujR) 2 , equation ( |3.2| ) 
reduces to a Bessel equation with the solutions 

0(r) = a.J l+2 (p)/p 2 + (3Y l+2 (p)/p 2 . (3.5) 

In the inner region, p < p <C 1, the last term in ( |3.4| ), Po/% 3 can be neglected in 
comparison with the typical^ "potential" term ~ 1/x 2 . Equation Q3.4| ) becomes 

d 2 <j> (l + x 2 )rf0 1(1 + 4) X 2 

+ ITT, ^ ( P = - ( 3 - 6 ) 



dx 2 x(l — x 2 ) dx 4x 2 (l — x 2 ) 4x(l — x 2 ) 2 

1 See 2.306 of jL7|. A complete absorption analysis based on the exact solution has been done in fli| . 

2 The term with I is an example of such a typical term. However, even for I = the "potential" 
characterizing the geometry is proportional to 1/x 2 - this can be seen by converting the equation into 
the Schrodinger form. 
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It has four regular singular points (at x = ±1, 0, oo). Unlike the case of a second order 
differential equation with three regular singular points which can always be cast into 
the form of hypergeometric equation by performing a suitable conformal transformation, 
such an equation has no common ancestor and should be studied on its own. The 
equation, however, contains a parameter A making it possible to study low- and high- 
temperature expansions even though the exact solution remains unaccessible. Finally, in 
the matching region, p <C (uR) 2 <C p <C 1, we have 

0(r) -> ap l /2 l+2 (l + 2)! - (32 l+ \l + 1)!/V +4 . (3.7) 

Anticipating that for ojR <C 1 (3 — > 0, the absorption probability can be written as 

P{1) = 27i\p 4 /\a\ 2 (3.8) 

To study solutions of ( |3.6|) , it will be convenient to eliminate the first derivative term by 
introducing u(x) = x~ x l 2 \Jx 2 — 10(x) and rewriting (|3.6| ) in the form 

u" = (\ 2 F{x) + G{x))u, (3.9) 



where 



1 . x 4 + Qx 2 -3 1(1 + 4) 
"77777 ^77-^ > G ( x ) = "777777 7777 + 



Ax{\ -x) 2 {l + x) 2 ' v ' ' 4x 2 (l-x 2 ) 2 Ax 2 {l-x 2 )' 



In the next subsection, we shall find solutions to equation (|3.9|) in the form of asymptotic 
series in 1/A. 

3.1 Low-temperature expansion I: Langer - Olver method 



We are interested in constructing a uniform approximation to the solution of ( |3.9| ) on the 
interval [0, 1], the endpoints of the interval being regular singular points of the differential 
equation. Our expansion will be built around x = 0, where the function F(x) has a simple 
poleQ. The idea behind the method of Langer and Olver is to introduce new variables 
in which the behavior of the solution close to the singularity would be fully described 
(for large A) by a known function, and treat the rest as a perturbation. In our case the 
situation is more complex since the interval of interest contains two singular points, and 
we have to prove that the functional coefficients of the expansion converge everywhere on 
[0, 1] - in particular, that they converge at x = 1. The convergence can be established by 
analyzing the singularities of the functions F(x) and G(x) in ( |3.9| ) fH |. We also confirm 
it by an explicit calculation. 

3 The expansion cannot be constructed starting from x = 1, since in that case the error-control 
function diverges (see remarks on p. 205 of 
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By making a conformal transformation ( = —p 2 (x), where 



1 11 + y/x 

p(x) = - arctan \fx H — log = (3.10) 

2 4 1 — Jx 



and introducing a new function W(x) = x 1//4 y / p(a;)/(a; 2 — i)u(x) = x 3 ^ 4 ^Jp(x)(p(x) we 
transform equation ( |3.6| ) into 



■/ ll fA 2 z/ 2 -l . y>(Q l 



^-\4c ^ + Tr- (3 - n) 

where i/ = / + 2, and 

^ W) = ^ + ^ 2 - 1 )(^ + ^ 2 - 1 ). (3.12, 
16jr(:r) 16a; 

The function is analytic at £ = (x = 0) and is treated as a perturbation. 

Then the formal series solution of ( |3.11| ) corresponding to the wave propagating from the 
infinity to the horizon is given by 

Wix) = { H &( X pW f A n(x)X- 2n - ^H®(\p(x)) ±B n (x)\-^ , 

(3.13) 

where Hjp(z) are Hankel functions. The coefficients A n , B n are independent of A and 
can be calculated successively using the following equations [0 

A n+1 (() = uB n (() - (B f n (() + J V(C)S»(C)dC, (3-14) 

B n (0 = -A' n (C) + ^ J W)Mt) -("+ 1/W)} ^ , (3-15) 

with A = 1. The formal series in ( |3.13| ) provides a genuine asymptotic expansion under 
the condition that the successive integrals in ( |3.14| ), ( |3.15| ) are convergent 0. As can 
be readily seen f\, this is indeed the case for our ip(x), and therefore Olver's theorem 
guarantees that (|3.13|) is an asymptotic expansion approximating the solution of (|3.11|) 
uniformly for x G [0,1]. In practice, it is far more convenient to use the following 



expressions for the coefficients A n (x), B n (x) (see Chapter 12, Exercise 5.2 in [14 



A n{X) = }^ h—\ A 2n-k{X) , (3.16) 

fc=0 P \ x ) 



4 Precise error bounds on partial sums of ( |3.13| ) can be found in Chapter 12, § 4 of j[4|). Conditions on 
ij) sufficient (but not necessary!) for the convergence of the error-control function are stated in Exercise 
4.2, Chapter 12 of [fj. 

5 For x — > 1 we have B n ~ 1/loge, A n ~ O(l), where e = \x — 1|. 
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2n+1 i k auiv\ 

B n {x) = i £ -0-J-A 2n+1 _ k {x) , (3.17) 

where 

/g^; 2 _|_ ^p, J 
— — 7= ^cfe, (3.18) 

and a fc (z/) = (4z/ 2 - l 2 )(4z/ 2 - 3 2 ) . . . (4z/ 2 - (2k - l) 2 ) /8 fc A;!. The normalization constant 
N^XJ) in ( TO , 

Ni(X, I) = e *Aiog2/2-i ff /4+iAV8 g ^(i) _ * li m p( x )5 n (x)A- 2 ") , (3.19) 



n=0 




is chosen in such a way that for x — > 1 we have 0(x) = x 3 ^ 4 W/ yp(x) (1 — x) 
Then in the limit i->0we obtain 

^UW^, (3 , 0) 
where 

oo 

iV (A, = £ A- 2n (A n (0) - 2uB n (Q)\~ 2 ) , (3.21) 

n=0 

which can be compared with aJi +2 {p)/ p 2 in the region p <C 1 to get the coefficient a. 
The absorption probability is calculated in a standard way, and we find 

where Po(l) is given by Q2.1|) . Calculating the first few coefficients, we get So(0) = 0, 
lim a; ^ip(x)So(a;) = (z/ 2 -l)/2, Ai(0) = (z/ 2 -l) 2 /8 (see also Appendix A). Then the ratio 
/z(A) = P{1)/Pq(1) is given by ( |2.3|) . Higher order terms can be calculated as a result of 
a lengthy but essentially straightforward procedure using the recurrence relations ([3.16 - 



3.18| ). The result ( |2.3| ) therefore potentially accounts for all power-like terms in the A ^> 1 



expansion. From our treatment of the Dl — D5-brane system in Section 3 it will become 
clear that the Langer - Olver expansion approximates the exact solution reasonably well 
even for A ~ 1. In fact, the first few terms of the asymptotic expansion ( |3.13|) coincide 
with the Frobenius series solution of ( |3 . 1 1|) at x — 0. 

The following assertions regarding the asymptotic expansion (|2.3|) appear to be true: 

(1) It is an expansion in powers of 1/A 4 . In particular, the 1/A 2 correction vanishes, in 
disagreement^ with g, 0, |K)|, 



6 In ||, the potential was approximated by its Taylor expansion near one of the singular points. Such 
an approximation is non- uniform and cannot give an accurate matching. The approximation used in Q , 
flCj is essentially of the form f(x) = x a ° (1 — x) ai , where ao,i are indices of the singularities at x — 0, 1, 
respectively. This approximation provides the correct result for P at A = oo (absorption probability at 
zero temperature) but it cannot be used to calculate finite-temperature corrections. We remark that the 
results of H , [flOl disagree with those obtained in M . 
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(2) For the s-wave [y = 2) power-like corrections in ( |2.3| ) vanish. 



Being an asymptotic expansion in powers of 1/A, formula ( |2.3D cannot account for the 
corrections of the form exp (—A). In the next section we show that such terms do indeed 



exist. Unfortunately, an efficient recursive scheme for their calculation is not known (19 
Our attempt to construct it seems to face considerable technical difficulties. 

3.2 Low-temperature expansion II: Reduction to a hypergeo- 
metric equation 

Since we are interested in studying solutions in the interval [0, 1], it seems plausible that 
singular points x = 0, x — 1 play the major role for such solutions, whereas the singular 
point x = — 1 can be treated perturbatively. Indeed, neglecting the term with G(x) for 
a moment, we can write ( pj.9| ) as 

U = ~Ax{l-x) 2 r 1 {x)' (3 ' 23) 

where 7](x) = (1 + x) 2 varies smoothly from 1 to 4 for x G [0, 1]. If we treat rj as a 
constant, equation ( |3.23|) can be solved exactly in terms of a hypergeometric function, and 
therefore one may hope that the solution of the original equation ( |3.9| ) can be successfully 
approximated by a series based on hypergeometric function. We shall now try to make 
these ideas precise. 

First, we introduce a new variable, ((x) e [0, 1], and make a conformal transformation, 

C = tanh 2 p(x), (3.24) 
p(x) = (arctani/x + arctanliA/x) /2 , (3.25) 



U = u^JdC/dx, (3.26) 
to write equation ( |3.9| ) in the form 

u " = - 4((i-cy u + vu > (3 ' 27) 

where derivatives are now taken with respect to ( and 

V = x 2 g(x(C))-l{x,Q ■ (3.28) 

Here x = dx/d( = y/x(l — x 2 ) cosh 2 p(x)/ tanhp(x), and 

r ^ x 3 fx\ 2 

{x,(} = 

x 2 \xj 
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denotes the Schwarz derivative. We have 

{ J^(l + x 2 /h + 0(x i )) forx^O, 
p i x \ — J (3 29) 

|-Jlog(l-x)+log2/2 + 7r/8 + 0(x-l) for x -»• 1, 

and so 

f x - 2x 2 /3 + 0(x 3 ) for x -»• 0, 
C x = 7 \ ; ' 3.30 

1 1 - 2v/T^e-"/ 4 for x -»• 1. 

Accordingly, the potential P(C) becomes 

f (z/ 2 - 1)/4C 2 + 0(1) for x -> 0, 
P(C) = (3.31) 
I -1/4(1 - C) 2 + 0(1) forx->l. 

Introducing the "regularized" potential, Vr = V — (y 2 — l)/4£ 2 + 1/4(1 — () 2 , we can 
write equation ( |3.27|) as 

u ''-{-W L W + ^-mT L o' + Vb ) u ' (3 ' 32) 

where Vr is a regular function for £ G [0, 1]. We shall look for solutions of the equation 
(13.321) in the form U = Uq(1 + h(Q), where Uq obeys 



4C 2 (1 - CW + (A 2 C + C 2 - (v 2 - 1)(1 - C) 2 ) ^o = 0, (3.33) 
and for the error-control function h(Q one can derive the following integral equation 



i C 



HO = J V R (1 + h(C))U 2 (0 J jj^dt . (3.34) 
The solution of ( |3.33|) is given by 



, . , . i-ix , . i+w ( v — iX , . v — iX , . \ 
Z7o(C) = (1 - C) — (-C)— aFi — s~ + 17+ v , — s~ + »7- " 5 1 - *A; 1 - C , (3-35) 



where tj±(i/) — ( 1 ± V^ 2 — lj /2. Returning to the x variable, the solution normalized 
to (1 — x)~* A / 4 for x — ^ 1 is 

O ( X ) = C A x 3/4 cosh iA p(x)tanh' /+1/2 p(x) 

2P1 (-iA/2 + z//2 + r? + , -iA/2 + i//2 + 1 - iX; sech 2 p(x)) , (3.36) 

where C\ = exp (iAlog2/2 — inX/8). For small x (recall that x = p^/p 2 ) one has 

C7 A r(z/) r(i - zA) 



p[, T (-iA/2 + i//2t7+(z/)) T (-iA/2 + u/2 + 7j_(i/)) 
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(3.37) 



This matches (for uR <C 1) the outer region solution, aJ u (p)/p 2 ~ ap l /2 l+2 Y{l + 3), with 

= 2^c A r(i/ + i)r(i/)r(i-»A) 

r(-iX/2 + p/2 + V+ (u))r(-iX/2 + p/2 + rj^u))- l ' ' 

The absorption probability ( |3.8|) can be written as 

2ttA |r {v/2 + r) + (u) + i\/2) T (i//2 + r?_(z/) + iA/2) | 2 /cur \ 2v 



\v\{y- l)!] 2 |r(l + zA)| 2 V 2 



(3.39) 



where we used |r(x + = \T(x — The zero-temperature limit (corresponding to 
A — > 00) gives (see Appendix ||) 

P = P (l + 4zV - 1)/3A 2 + ■ • •) , (3.40) 

where Pq is given by (|2.1|) . We must not forget, however, that there are corrections to the 
result ( p.39| ) coming from the function h(() in ( |3.34| ) . Standard analysis of the Volterra 



equation Q3.34D suggests that h(() is indeed a small correction (properties of Vr and U 
ensuring convergence of corresponding integrals at the endpoints) of the order of 1/A 2 . 
Therefore only the zeroth order term in the expansion ( |3.40| ) can be trusted and thus 
there is no contradiction with our main result Q2.3p . The significance of ( ^.39[ ) is twofold: 



it shows that the absorption probability must contain contributions of the form exp (—A), 
and that it has a functional form expected in a field theory (compare (|4.4| )). We can try 
to guess the exact expression for the absorption probability by choosing the ansatz of the 
formF] 

_ 2ttA |r (xi + iA/2) r (X2 + iA/2) | 2 (cur \ 2 » 

H(i/-1)!] 2 |r(l + zA)| 2 \ 2 ) [ ' ' 

Comparison with the asymptotic expansion ( p.3| ) to the order 1/A 2 gives 



xi, 2 = \ v+l± ^(l-^)/ 3 j /2 . 

However, the coefficient in front of 1/A 4 term turns out to be ?>2v(i> 2 — l)(4z/ 2 — 1)/45 which 
does not match ( [2.3|) . This suggests that the naive thermal ansatz ( p.41|) is incorrect. 

3.3 High-temperature expansion 

Two independent solutions of ( |3.6|) are 

<f>±{x) = (l-x) a± F(x), (3.42) 



7 For the Dl — D5 system this method gives the exact thermal factor in (4.4) immediately. 
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where a± = ±z'A/4 are the roots of the indicial equation and F(x) obeys 
d 2 F x 2 {l + 2a) + 2ax + ldF 1(1 + 4) 



2\ 



dx 2 x(l — x 2 ) dx 4x 2 (l — x 

X (l +X ) a(ff _l) +(T (1 +X 2 ) A 2 

x{l-x) 2 (l+x) 4x(l-x 2 ) 2 1 ' 

The solution corresponding to the incoming wave at the horizon is 

4>{r) = A(l - r 2 /r 2 )- iX ^F(r) , (3.44) 

where A is an arbitrary constant and F[r) is supposed to be regular at r = Tq (x — 1). For 
A 1 the approximation scheme can be constructed applying a standard perturbation 
theory to the equation (|3.43|) . In the zeroth-order approximation Q3.43j ) reads 



x 2 (l - x 2 )F" - x(l + x 2 )F' - /(/ + 4)/4 F = , (3.45) 

and has solutions 

F 1 (x) = x 2+l/2 2 F 1 (l + 1/4, 1 + 1/4; 2 + 1/2; x 2 ) , 

F 2 (x) = x' l/2 2 F 1 (-l/4,-l/4;-l/2;x 2 ) . (3.46) 

We have to distinguish between even and odd I. 
Odd /: 

The solution regular at x — 1 (i.e. the one without logarithmic singularities) is given by 
F(x) = Fi(x) + QF 2 (x), where 



r(2 + |)r 2 H; 



r(-|)r 2 (i + |) 



Matching this with the outer region's solution gives a = 2 l+2 [l + 2)!Cj/ p , and the 
absorption probability 

where Pq(1) is given by (|2.1|). 
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Even I: 



In this case the solution regular at x — 1 is Fi(x). For / = 4n, n — 0, 1, 2, it is given 
by 



(2n)\ n 



where P„(z) are the Legendre polynomials, and we have -^(x) 



-i/2 for 



X 



0. For 



Z = 4n + 2, n = 0, 1, 2, ... the solution can be written in terms of the Meijer function |20 

F 2 (x) =G™ (x 2 ' , ' 1 l ) . 

2K J 20 \ -n- 1/2 n + 3/2 • - J 

x -i/2 /r 2 (l + 1/4) for x — > 0. Accordingly, the absorption probability 



We have F 2 (x) 
is 



p = p o(/) C^y- (!±i)!! = m (I±D1 



for I = 0,4,8, ... and 



p m 2 '+ 3 {i + iy 2 v\i + i/4) = p (2+r\ a + 3 (/ + i)! 2 r 4 (i + //4) 



>(0 



|2p4 

7T 



2Z+3 



(0 



V u 



■) 



7T 



(3.48) 



(3.49) 



for I = 2,6,10.... Thus for the s-wave the absorption cross section cr = 32ir 2 P (0) / u 5 
can be written as o"o = A(ro, R)f(X), where A(tq, R) = n^r^R 2 is the horizon area of the 
three-brane^ in the near-extremal limit, and /(A) interpolates between 1 at A = and 
7rA 3 /8 at A = oo. 

Corrections to the lowest-order results (|3.47|) , (|3.48j) , (|3.49|) can be obtained by applying 
regular perturbation theory to the equation (|3.43|) . For example, for the s-wave we have 



F(x) = 1-— lot 



l + x\ A 



32 



l0£ 



1 + X 



+ 8 log 



l + x 



4Li 2 



1 — x 



Correspondingly, the absorption cross section becomes 
a = n 6 R 8 T 3 



(tt 2 + 12(2 -log 2) log 2) iv 2 



+0(A 3 ) . 
(3.50) 

(3.51) 



48tt 2 T 2 

The frequency- dependent term in (|3.51| ) can be interpreted as a correction to a hydro- 



dynamic description of the strongly coupled Yang-Mills plasma p2 



'This is reminiscent of the universal result for the black holes [Ell 
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4 Absorption by the Dl / D5-hra.nes 



The bound state of Dl and D5 branes has been intensively studied as a string theory 
model for a five- dimensional black hole, and later in the context of AdS/CFT corre- 
spondence (discussion and references can be found in IJ). The non-extremal absorp- 
tion probability contains a thermal factor which can be calculated exactly using a (su- 
per) gravity description in the "dilute gas" regime. Here we would like to calculate this 
temperature-dependent coefficient by using L anger- Olver's approximation. We shall see 
that power-like terms at any given order of the low-temperature expansion are indeed 
correctly reproduced, but the corrections exponentially suppressed at low temperatures 
are not captured. We start with the "near region" equation (3.7) of M, where we shall 
put Tl = Tr = Th = T for simplicity and use notations x = r 2 /r 2 , X = u/2nT: 

_ ^ + 2) + A 2 = Q ^ 

dx 2 (1 — x) dx 4x 2 (l — x) 4x(l — x) 2 

Introducing a new function u by u = y/x — 1 (j) one can write (|4.1|) in the form 

u" = (\ 2 F + G)u, (4.2) 



where F = -l/4x(l -x) 2 , G = 1(1 + 2)/4x 2 (l - x) - 1/4(1 -x) 2 . Equation Q has an 
exact solution given (up to an arbitrary constant) by 

= X ' l/2 (1 - x)~ iX/2 2 F 1 (-(I + iA)/2, -(I + zA)/2; 1 - zA; 1 - x^j (4.3) 

which leads to the absorption probability [EJ 



{ ) l\ 2 (l + l)\ 2 \ 2 J 



p2 ( 1+2 iX 
\ 2 2 



r(i-iA) 



(4.4) 



The plane wave absorption cross-section is then obtained by multiplying ( f4.4| ) by Att/uj 3 . 
In the low temperature limit (A — > oo) we have 

v(y 2 -\) v(y 2 -\)^-Vv 2 -^ + 1\) 
P(/)/P (/) - 1 + + — 

t/(t/ 2 - l)(i/ - 2){y - 3)(35z/ 4 - 14z/ 3 - 80z/ 2 + 314z/ + 465) 
+ 5670A 6 

+ 0(l/A 6 )+0(e~ A ) , (4.5) 



where v = I + 1, the absorption probability at zero temperature is given by 

P (1) 4?f2 ( 21+2 = K 2 uj 4 r 2 r 2 (uj 2 ri r 5 ) 21 

oU Z! 2 (/ + l)! 2 l 4 4 2 4 'Z! 2 (/ + 1)! 2 ' 1 J 
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and we have used T = r^jl-nrxr^. Notice that in addition to the power-like corrections 
coming from gamma-functions in the numerator of (FOj), expression ( |4.5| ) also contains 
corrections of the form exp (—A) produced by \T(l—i\) | in the denominator (see Appendix 
|B|). In fact, for I = there are no power-like corrections at any order: 



p(0) = \2^J T = Po(0) ( x + 2e ~ 



-A 



(4.7) 



4.1 Low-temperature asymptotic expansion 

Using Langer-Olver's method, we would like to construct an approximate solution of 



(4.1) in the form of asymptotic series in powers of 1/A 2 , and compare it to the large A 
expansion of the exact solution. We shall use the notations of Section [| throughout, with 
the exception of A = uj/2ttT and v = I + 1. Applying the conformal transformation 



we can write (14.21) as 



where the function 



W 

tp(x) = 



( = — arctanh 2 \/^ , 

arctanh 1/ ' 2 v / x 

w = — rn—, u 
x l / A yfx=l 



1 -4z/ 2 (1 - x)(l - 4z/ 2 - x) 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



16C(x) ' 16a; 

is analytic at x — 0. The formal series solution of ( f4.10|) representing the incoming wave 
is given by 

W(x) 



/ttA 



. ./— arctanhy/q? I ffff, (A arctanh\/x) A n (x)X 2 " 
Ni[X, I) V 2 ^ 



n=0 



A 



—=H^\(\ arctanhv^) £ B n {x)\~ 2n 

arctann^/ x n=o 



(4.12) 



where the normalization constant Ni(X,l), 



n=0 



N\(X, I) = e Alog2 - k/4 Y AJ1) + - lim arctanh v ^5 ri (x)A 



-2d 



is chosen in such a way that 



x 



1/4 arctanh- 1/2 v^ W -> (1 - x) 



-i\/2 



(4.13) 



(4.14) 
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for x — > 1, and thus can be identified with the exact solution ( [4.3| ). Since for x — > 

- A" 2' i' Z! 



A^A^ 1 



;r 



-Z/2 



where 



iV (A,/) = X! A" 2 " (X(0) - 2z/ J B„(0)A- 



n=0 



(4.15) 



(4.16) 



one can match with the outer region's solution aJi + i(p)/p in the overlapping region 
p <C 1, and obtain the absorption probability 



P(Z) 



4vr 2 



/ Au;r N 


2Z+2 

) 




2 = Po(0 


iVi 


V 4 , 




N 





(4.17) 



The recursion relations for the coefficients ( |3.16| ) - ( |3.17| ) are now 

Au 2 - 1 + x T 



A fc+1 (x) =iy/x(x- l)d x A k + 



16x\/x 



One finds 



S (x) 



4z/ 2 



Sarctanh^^/x Sarctanh-y/a; 

so that -Bq(O) = Z(Z + 2)/6, and lim^i arctanhy / iri?o(a;) = z^ 2 /2 — 1/4. Limiting values of 
a few higher-order coefficients can be found in Appendix |A|. Substitution into (|4.17|) gives 
precisely (|4.5| ) as expected. Due to the simplicity of the function ip(x), the calculation 
of coefficients A n , B n is straightforward, and therefore our result can in principle be 
extended to all orders in 1/A 2 . However, the asymptotic expansion Q4.12|) cannot capture 
exp (—A) correction terms which in the s-wave case would lead to an erroneous conclusion 
that -P(O) = -Po(O) for non-zero temperature. 

In Figure ([[]) we plot the exact and the approximate solutions of (|4.1|). Figure (g) 
shows the ratio P(1)/Pq(1) (exact and approximate results) as a function of 1/A. It is 
clear that the Langer-Olver approximation which is supposed to work well for large A 
is in fact quite good even for A ~ 1. We remark that the Pade- improvement of the 
low-temperature asymptotic series ( |4.5|) does not seem to be particularly helpful in this 
case. 



1 + 



l-4z/ 2 



(4.18) 



(4.19) 



4.2 High-temperature expansion 

The leading term in the high-temperature expansion (A 1) of ( |4.4|) is 

'4irT H \ 2l+1 T 4 (l + Z/2) 



71 



(4.20) 



16 



For the s-wave the absorption cross section is given by uq = 2ir 2 rorir5f(X), where /(A) = 
(jrX/2) coth (ttX/2) interpolates between 1 at A = and ttX/2 at A = oo. Corrections can 
be readily obtained by solving equation ( flip perturbatively in A < 1. For the s-wave 
we get 

4>(x) = (1 - x)- lX ' 2 (l - ^Li 2 (l - x) + ■ ■ ^ . (4.21) 
Correspondingly, we have for / = 0: 

P(0)/P o (0) = 2/ttA + ttA/6 + 0(A 2 ) . (4.22) 

For I = 1 the result is 

PiX)IW) = + 1 i 1 ~ Ya) + 0(A) • (4 ' 23) 



Of course, formulae Q4.22 ), ([4.23|) coincide with the A C 1 Taylor expansion of the exact 
result (OI). 



5 Field theory calculation 

In field theory, the absorption cross section at finite temperature can be calculated using 



the method of thermal Green's functions (see e.g. [23|). The absorption cross section is 
determined by a discontinuity of the two-point function of the operator 0\ to which the 
Zth partial wave of the dilaton field couples, and is given by the zero-momentum limit of 
the spectral density, 

of = (5.1) 

where pf(u, k) at finite temperature in the imaginary time formalism is proportional to 
the discontinuity of the (analytically continued) Euclidean thermal two point function 
n(r,r), 

p f (^ k) = jj d 3 xdte~ ik * +iuJt [n(r, t-ie)- fl(r, t + ze)] . (5.2) 

Here the function 5(r, t) denotes the analytic continuation of IT(r, r) to the real time 
t = —ir, T denotes the flux, and Q is a coefficient which depends on the parameters of 
the system under consideration. 

5.1 D1/D5 system 

The absorption cross section is given by ( |5.1| ), (|5.2|) , where Tl(r,t) is the (analytically 
continued) two-dimensional thermal Green's function, 

nr t)~ c ° ( nn x hL ( nTR y hR 

) i2h h +2h R Ui n h7rT L (t + x) ) \smh7rT R (t - x) J ' 1 j 
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Figure 1: D1/D5 system: exact solution Re 4>(x) (solid line) and a partial sum (with all 
A s , B s equal to zero for s > 0) of the asymptotic series for A = 1 (dashed line). 




Figure 2: Dl/DB system: ratio P(l)/P (l) versus 1/A for I = 1: exact (solid line), 
high-temperature approximation ( |4.23| ) (dashed line), sum of the first three terms of the 
low-temperature asymptotic expansion ( |4.5| ) (dotted line). 
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The integral in ( |5.2| ) factorizes nicely by introducing x± — t ± x, and can be computed 
exactly f24" |. This factorization is a property of a two-dimensional theory and does not 
hold in higher dimension. To illustrate the approach we shall use in the d = 4 case we 
compute the D1/D5 absorption cross section <Jq by integrating over t and x directly 
without factorization. One can also show that the low-temperature asymptotic series 
( }4.5|) can be obtained by expanding Tl(r,t) in powers of T. The spectral density for the 
s-wave is given by ( |5.2|) with Q = 2Txr\r\, T = to and 



ti{r,t) 



ir*T 2 T 2 



sinh 7rT^(t + x) sinh TcT R (t 



x) 



Integrating over time we get 



dxe 



-ikx 



e -iu>x T 2 ( 27 riT R coth (2ttT r x) 



to 



sinh 2itT r x 



e iux T 2 (2-niT L coth {2nT L x) + u) 
sinh 2-kTlx 

Evaluating residues, we obtain the spectral density 

n 2 (cu 2 - k 2 



Po (w, k) 



, (u + k) , (u - k) 
coth v ,„ ; + coth v 



4T, 



ATr 



The absorption cross section is then obtained by taking the limit k 

- 1 



in (|5~ 



a = a f(uj,T L ,T d 



R 



12 2 



(e^L - l) {e^R - l) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



where Oq is the cross-section at zero temperature. Similar calculations can be done for 
higher partial waves. For all of them one observes a perfect agreement with the results 
obtained from gravity ||, [§, @, @. 



5.2 Non-extremal three-branes 

To the leading order in 't Hooft coupling, the terms contributing to (Oi(x)Oi(O)) are a 
bosonic term IT^ 08 (for all / > 0), a two-fermion term U 2e (for / > 1), and a four-fermion 
term Ilf (for / > 2) identified in [p5| : 



n 



T=0, 

£,bos ' 



X 



^C+i^Tl^^ '^ (^ A W^ A W + ^ 2 A(x)d 2 A(x) ) (5i 



where A(x) = l/47r 2 \x\ 2 is the free scalar propagator. The term proportional to d 2 A(x) 
does not contribute to the discontinuity and will be omitted. At finite temperature the 
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bosonic contribution to the two-point function will be given by ( |5.8| ) with A(x) replaced 
by its thermal cousin. Borrowing the expression for the thermal propagator in position 
space from P5| , 

„ , N T sinh27rrT 
Gb\t, t) 



Anr (cosh 2irrT — cos ItxtT) 



where r = we obtain 



rp . . k 1 uj 21 N 1+2 (1 + 2)U + 3) , x 

n ^ (r ' r) = 3,^11(1 + iy. M(r ' r) 



where 



M(r, r) = ^ [(<9 2 r £ B ) 2 + («9 2 r £ B ) 2 + 2{dl T G B f + 2(d r Q B ) 2 /r 2 
For the s-wave we obtain (see Appendix |C| for details) 



2\2 



7lk(tU 2 - k 2 ) 



32 



2T 
1 + — log 



1 — e 2T 
1 — e 2T 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



Taking the limit k — > we get 



<7 n = coth — — = cr n coth — — , 

o 32n AT AT 1 



(5.13) 



where n 2 N 2 = 4ir 5 R 8 . 

For partial waves with I > in addition to (|5.2| ) one has to take into account contributions 
of the two- and four-fermion operators. The contribution of the two-fermion operator is 



nT ^ 2i iV' +2 (/ + 2)(/ + 3) 



A^ti,2e(r,r) 



3ttW(Z- 1)!(Z + 1)! 
where 

M lj2e (r, r) = g l B l [d 2 „Q B (d r g F y + d 2 TT g B (d T g F ) 2 + 2d 2 T g F d r g F o T g F 

and the free fermionic propagator at finite temperature is given by 

sinh irrT cos 7tTt 



Gf(t,t) 



T 



2nr (cosh 2-nrT — cos 2htT) 



(5.14) 



(5.15) 



(5.16) 



Here we only consider the I = 1 case. For general I the calculation is similar but more 
tedious since the four-fermion operators also contribute. Expanding Aii^eir, T ) in powers 
of T and evaluating the discontinuity we get 



a 



1, field th 



6144vr 7 /2 



'ttT 



V to J 



(5.17) 
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Comparing this result with the one obtained from gravity. 



<7i 



K 3 N 3 UJ 7 

6144vr 7 /2 



l + 4£ 



10 



O (e-^ T ) 



(5.18) 



we observe that crf floldth contains terms proportional to T 2 and T 6 which are absent in 



6 Black £>-branes and non-extremal M-branes 

In this Section we demonstrate that the M- and black p-brane backgrounds can be treated 
by the approximation methods used in Section 2 provided the correct identification of 
the expansion parameter is made. 



6.1 Non- extremal M-branes 

Non-extremal M-brane solutions were found in [[27[ . In the notations of [28 the M2-brane 
metric is given by 

ds 2 n = H- 2 l'\r) (-f(r)dt 2 + dxl + dx 2 ) + H^r) (J^ + r 2 dQ 2 ^ , (6.1) 

where H(r) = 1 + i? 6 /r 6 , f(r) = 1 — r^/r 6 . The radial part of the Laplace equation for 
a minimally coupled massless scalar propagating in the background ( |6.1|) becomes 



d 2 f (5x 3 + l) df 1(1 + 6) X 2 ( p% \ 

dx 2 2z(l - x 3 ) dx 4x 2 (l - x 3 ) J (1 - x 3 ) 2 V {ujR) 6 x 3 J J ' 1 ' ' 

where x = r^/r 2 , A 2 = (uR) 6 /Ap^, p = ur. In the "inner" region p <C 1 equation (|6.2|) is 

„, , / A 2 1(1 + 6) 5-46x 3 + 5x 6 \ . . , , 

w {x) = [-JT^y + 4x 2 (i-x 3 ) + i6x 2 (i-x 3 ) 2 ) w{x) (6 ' 3) 

where w(x) = f(x)x~ 1 ^y/x 3 — 1. Asymptotics of ( |6.3| ) can be constructed using Langer- 
Olver's method. Here, however, the parameter of the expansion should be redefined 
according to A 2 — > A 2 = A 2 - 31(1 + 6)/4 to ensure the convergence of the error-control 
function^ 0. 

9 This redefinition is similar to the one performed in quantum mechanics where in quasi-classical 
approach to the central-force potentials a substitution 1(1 + 1) — > (1+ 1/2) 2 is made "by hand" to achieve 
the correct asymptotic behavior of the solution. 
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In the case of M5-branes the metric is given by 

< = zm^±M±J^±M±M±M + ^/ 3(r) + ^ ? (fU) 

where if(r) = 1 + R 3 /r 3 , f(r) = 1 — r^/r 3 . An analogue of the equation (|6.2j ) is 

tl ^~ 2 ) d f ^ / + 3 ) / | A2 fi | ^ -V Q, (6.5) 
dx 2 x(l — x 3 ) dx x 2 (l — x 3 ) x(l — x 3 ) 2 \ (cuR) 3 x 3 J 

where x = r /r, X 2 = (ujR) 3 / p . In the "inner" region, i.e. for p < 1 equation ( |6.5| ) 
becomes 

,n \ ( ^ 1(1 + 3) x{8-x 3 )\ , . 

w{x) = {-xjr^xW 2 + x^k + 4(1^ J w ^ (6 - 6) 

where w(x) = f(x)x/(x 3 — l) 1 / 6 . As before, an asymptotic expansion with the convergent 
coefficients can be built, provided the redefinition A 2 — > A 2 = A 2 — 4 — 31(1 + 3) is made. 



6.2 Black p-branes: 

The metric is given by 

—f(r)dt 2 + dx ■ dx 



ds 2 l0 



H P (r) 



+ \HJr) 



7F) +r2dnlp ] 



(6.7) 



where H p (r) = 1 + R 7 p /r 7 p , f(r) = 1 — r p /r 7 p . Let p = cur. Then the radial part 
of the Laplace equation for a minimally coupled massless scalar can be written as 



dp 2 



a + 1 + 



a(uoR) a (4: - a) 
(p a -p a ) ' 2(p- + (uR) a ) 



+ 



ldf p a - 2 l(l + a) p a (p a + (ioR) a ) 

~J+ I „ _n\9 J " 



p dp 



P - Po 



(P a ~ Pof 



(6. 



where a = 7 — p, p = ur . Solutions in the "outer" region (characterized by p <C p , 
p (ujR) a ^ a ~ 2 ^) are given by 



/(r) = ap- a / 2 J l+a/2 {p) +(3p- a / 2 J^ a/2 {p) 



(6.9) 
(6.10) 



for odd a and by 

f(r) = ap- a J l+a/2 (p)+(3p~ a Y l+a/2 (p) 

for even a. To study solutions in the "inner" region, we introduce^ x = r 2 /r 2 for p- 
branes with odd p and x = r /r for p-branes with even p. Then for the odd p-branes the 

10 The difference in treatment of even and odd p reflects our desire to avoid fractional powers of x in 
the potential. 
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equation ( |6.8| ) becomes 



dx 2 



+ 



ax 



a/2 



+ 



_ a(A-a)x a ' 2 (ujR) a /p a 
~2x 2x(l - x a l 2 ) ~ Ax(l + x a / 2 (ujR) a /p$) 
A 2 (. p a n \ „ 1(1 + a) 



dx 



x 3-a/2(l - x a/2y 



1 + 



(ujR) a x a / 2 



f 



Ax 2 {l-x a l 2 ) 



f = 



(6.11) 



where a = 7 — p and A = (uR) a /ApQ . For even p-branes we get 



dx 2 



+ 



a ax 
+ 



a-l 



a(4 - a)(uR) a x a /f% 



d£ 
dx 



+ 



4-a/ 



x (1 - x a ) 2x(l + (ujR) a x a /p%) 

— (i + pg/M)^) / - i! +a L / = 



x 2 (l — x a ) ' 



(6.12) 



where a = 7 — p and A 2 = (uR) a /p^ 2 . In the "inner" region, p C 1, equations (|6.11|) 
and ( |6.12| ), correspondingly, become 



d 2 f (a 2 - 6a + 4 - (a - 2) V 1 / 2 d/ 



A 2 



dx 2 



4x(l - x a / 2 ) 



-J- 



1(1 + a) 



rf 2 / a 2 - 6a + 2 - (a 2 - 8a + 2)x a a 1 / 



dx x 3 ~ a / 2 ( 1 - x a l 2 ) 2 4a; 2 ( 1 - ) 
A 2 /(/ + a 



/ = 0. (6.13) 



-J 



~J = 0. (6.14) 



rfx 2 2x(l-x a ) dx x 4 - a (l-x a ) 2 x 2 (l-x a )' 

Equations Q6.13 ) and Q6.14 ) have regular singular points at x — 0, 1, oo and a number of 
others on the unit circle. These equations are in this sense similar to the black three-brane 
equation and therefore can be treated by the approximation methods used in Section |3|. 



7 Conclusions 

In this paper we develop a reliable approximation scheme to compute high- and low- 
temperature corrections to the absorption cross section of a minimally coupled massless 
scalar in the background of near-extremal black branes. The absorption cross section for 
the non-extremal three-branes is explicitly computed using this scheme, and compared 
with the field-theoretic calculation. Our main results are summarized in Section 0. 

A number of conceptual and technical points remain to be clarified. To construct 
a well-defined scheme based on a hypergeometric approximation (Section [3. 2D or any 
modification thereof would be an important development by itself. It would allow one to 
control exponentially small corrections and therefore quantitatively compare field theory 
and gravity results ( |2.4|) - Q2.5Q for the s-wave at low temperature. For I — 1, it would 
be interesting to determine the origin of the disagreement between (|2.6| ) and (|2.7|) for 
T ^ u. In field theory, the spectral densities of correlators with I > can in principle 
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be computed exactly. In gravity, corrections to the leading high-temperature results for 
higher partial waves Q3.47| ), ( |3.48| ), ( |3.49| ) can also be determined. 

It may also be interesting to repeat our calculations for non-minimally coupled fields. 

The approach outlined in Section |] could be useful for studying the world-volume 
theories on M-branes, as well as lower-dimensional QCD (see Section 6.2 of [p]] and 
references therein). 

Finally, the results obtained in this paper may prove useful in determining transport 
coefficients of the strongly coupled finite-temperature gauge theory p2| . 



A Coefficients A n , B n 

Here we provide explicit expressions for the boundary values of the first few coefficients 
A n , B n defined by ( CT) , ( gT7[) . 
D1-D5 brane system: 

A = 1, 

Ai(0) = (32z/ 3 -12z/ 2 -32z/ + 3)/192, 

Ai(l) = -(16z/ 4 -8z/ 2 + 2)/128, 

A 2 (0) = (10240z/ 6 - 44544z/ 5 - 6800z/ 4 + 132480z/ 3 + 17800z/ 2 - 87936i/ - 4185)/737280 , 

A 2 {1) = (128z/ + 2176z/-4496i/ 4 + 3704z/ 2 - 807)/49152. (A.l) 

Note that B n below is defined as lim x ^ip(x)B n (x). 

B (0) = (^ 2 -l)/6, 

Bo(l) = (2^ 2 -l)/4, 

5i(0) = -(116z/ 4 - 345z/ 2 + 229)/1920, 

#i(l) = -(16z/ 6 + 64z/ 4 -134z/ 2 + 33)/768. (A.2) 
Black three-brane system: 

Aq = 1, 

Ai(0) = 0, 

= -(^ 2 - l) 2 /8, 

A 2 (0) = (795 + 32768z/-420z/ 2 -40960z/ 3 - 11760z/ 4 + 8192z/ 5 + 2880z/ 6 )/40960, 

A 2 (l) = (6801 - 39404z/ 2 + 9648z/ 4 - 2880z/ 6 + 320z/ 8 ) /12280 , 

A 3 (0) = 0, 

A 3 (l) = 1 146474829 - 1016514774z/ 2 + 1282758581z/ 4 
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- 553717340z/ 6 + 95818800z/ 8 - 5888960z/ 10 -24640z/ 12 J/1135411200. (A.3) 

B„(0) = 0, 

B (l) = (V-l)/2, 

Bi(0) = (z/ 4 - 5z/ 2 + 4)/5 , 

Bi(l) = (1344- 9152z/ 2 + 3520z/ 4 - 320z/ 6 )/15360, 

£ 2 (0) = 0, 

H 2 (l) = ^76248624- 46481 1632z/ 2 + 280281408z/ 4 

- 60890368z/ 6 + 3942400z/ 8 + 78848z/ 10 ^J /302776320 . (A.4) 

Higher coefficients can be rather easily obtained using Mathematica™ but the expression 
is complicated and will not be presented here. 



B Asymptotic expansion of \F(x + iy)\ for \y\ — » oo 

Using the Stirling expansion one can obtain the following formula: 

\r(x + iy)\ = y/tore-*W 2 \y\ x -V 2 {l + x{x - l)(2x - l)/12|y| 2 

+ a;(20a; 5 - 132a; 4 + 245a; 3 - 150a; 2 + 5x + 12)/1440|y| 4 

+ x(280x 8 - 4284a; 7 + 23046a; 6 - 50337a; 5 + 42735a; 4 

- 2331a; 3 - 10801a; 2 + 252a; + 1440)/362880M 6 + O(l/|y| 8 )} (B.l) 

Notice that for x — 0, 1/2, 1 the expression in the curly brackets is identically equal to 1 
(no power-like corrections). This is not a surprise since we have |T(ix)| 2 = n/x sinh nx, 
|T(l/2 + ix)\ 2 = n I cosli7ra; and |T(1 + ix)\ 2 = ttx/ sinli7ra;. Another useful fact is the 
product expansion of F(x + iy): 

r(x + iy)=^, (B.2) 

where 

oo r 2 ^-i/2 

t = rw n{ 1 + (^} ■ (R3) 

1 = y{ -IE + £ { \ - - tan" 1 } . (B.4) 



n=l 



n y x + n — 1 



25 



C Integrals required for a thermal field theory cal- 
culation of the absorption cross section 

Our approach to evaluating relevant thermal field theory integrals and their low- and high- 
temperature asymptotics can be illustrated by the following example. We are interested 
in computing the integral 



I(b,c) 



cosbx , / nb\ ixb 

— — — dx = - log 2 cosh — = - — 
x sinn cx \ 2c I 2c 



- log f 1 + e- nb/c 



(C.l) 



where the right-hand side is defined in the sense of analytic continuation. To calculate 
flC.lj ) exactly one can either take the limit a — > of the expression (2.5.47.1 (p. 469) of 



x 



' ' 1 cos bx 



x sinh cx 



dx 



r(q) 

(2c)" 




+ c 




(C.2) 



or evaluate the sum of residues with an appropriate prescription for the pole at x — 0. 
To obtain a "low-temperature" asymptotic, c — > 0, we can expand sinhcrc in series, 
and then integrate using residues. Taking c — > pushes all poles of sinhcx except 
x = to infinity, and thus our integration picks only the "polynomial" part of the exact 
result, —Tib /2c. Exponentially suppressed terms are not captured. This, however, is 
exactly what we need to compare the field-theoretic absorption cross-section with the 
one obtained from gravity in the form of a low-temperature expansion. 
The other limit, c — > oo, is more subtle. In this limit, the poles of sinhcx merge on the 
imaginary axis and one cannot approximate sinh cx by the exponent. One way to proceed 
is to rescale z = cx, expand cosbz/c for large c, and then integrate. This produces (in 
the zeroth-order approximation) the right answer, —log 2, assuming the regularization 
provided by the integral 

oo , 

x 2 a — 1 

- dx = 9Q _i r («)C(«) • 



sinh x 



The spectral density of the two-point function ( 5.10|) for I = is given by 



Po (w, k) 



K 



TT 2 L0k 



(P(2ttT, u — k) — P(2vrT, u + k)) 



(C.3) 



where 



P(a, b) = A(a, b) + B(a, b) + I(a, b) + J(a, b) + S(a, b) + T(a, b) , 
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a, t\ 3a 2 , . 7 , a 4 / cj 2 \ , . 3a 4 , . 7 . 

A(a,b) = —A i:2 (a,b) + —{l + — \A 2i2 (a,b) + —A 2A (a,b), 

,\ 3a 2 w „ , a 4 u . 3a 4 cj . 

5(0,6) = -—53,2(0,6) + — J B 1>2 (a,6) + — ^,4(0,6), 

r/ ,x 3a r acu 2 7X 3a 3 . 7 . 

AM) = - — 4i(°> & ) + -g-^3,3(a,6) 

t/ n 3aw 7X 3a 3 cj 7X 

</(a,6) = -g-^uM) + — -MM) , 

S(a,6) = ~ Y^ 3 ' 

T(a,6) = -^~ T 4, 

and the integrals involved are defined by 

. . IN r cos br dr r sin br dr 

Am ^ h) = /r^nl?V' B ^(a,b) = J- 



r m sinh™ ar 

cosh ar cos 6r ar . , . f cosh ar sin br dr 



f cosh ar cos br dr . f 

WM) = J rmsinh n ar , 6) = J 



r m sinh™ ar 



. sin6rdr fcosbrdr 



Calculation gives 



r . 7X tt6 (8 a 4 + 20 a 2 6 2 - 3 6 4 ) 2a 4 p 5 

VM) = ^77T + ^r~ 



4i(a,6) = 



360 a 7T 4 
6(-2a 2 + 6 2 )vr 2a 2 p 3 



6 a 7r 



2 ' 



tt6(6 4 - 8a 4 ) 12a 2 p 5 6a6p 4 , & 2 P3 
/3 < 3M) = 1200^ + ^^ + ^ + T 2 "' 



^1,3(0, 6) 



7r6 3 2p 3 26p 2 & 2 Pi 
6a 3 7r 2 7ra a 2 



r , , . 7r6 3 (20a 2 + 36 2 ) 2p 5 26p 4 ^ / 2 6 2 \ 
Jl < 5M) = 360a^ + ^ + ^ + S^ + ^r 3 



/ 6 3 26 \ / 6 4 6 2 \ 

+ TT— - + P2+ T7TT + ^ Pi, 
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T f h\ 7Fa i ^ 1 6a ^ 4 , ^3 . b 2 p 2 

15 24cr tt 6 7i z na 
T / ,\ ^a? -nab 2 nb 4 2a 3 p 4 4S 

, , L , vr6(24a 4 + 20a 2 6 2 + 36 4 ) 8a 3 p 5 2a 2 6p 4 ,~, r , 
A4 ' 2(G ' 6) = 360^ r? (C - 15) 



7rfe(2a 2 + fr 2 ) 4ap 3 , 2fep 2 

OCT 7T Z 7T 



„ , , , 7r6(88a 4 + 40a 2 6 2 + 36 4 ) 8ap 5 Qbp 4 ( 8a 3b 2 \ 

M > M = 360^ ^ ~ — -(3^ + ^J P3 

( b 3 46 \ 

+ 7HP2, (C.17) 



3a 2 7T 37T 



oa 7r a 



„ , 7N vr(88a 4 + 120a 2 6 2 + 156 4 ) 2p 4 26p 3 /2 6 2 \ 
Sl ' 4(a ' 6) = 360^ 73--^- 2 (^ + ^J P2 

R , 7r(8a 4 + 20a 2 6 2 + 56 4 ) 6a 2 p 4 2abp 3 

B 3t2 {a,b) = —— — , C.20 

120a' 1 71 tt z 

where p n (a,b) = Li n [exp (— irb/a)] is a polylogarithm jpCfl . We also have S3 = — vr6 2 /2, 
S 5 = vr6 4 /24, T 5 = vr6 3 /6. 
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